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MODIFIED NASA-LEWIS CHEMICAL EQUILIBRIUM CODE 
FOR MHD APPLICATIONS 

by 

R. A. Sacks, H. K. Geyer, S. J. Grammel, 
and E. D. Doss 

ABSTRACT 

A substantially modified version of the NASA-Lewis Chemical Equilibrium 
Code (1) has recently been developed. The modifications were designed to 
extend the power and convenience of the Code as a tool for performing com- 
bustor analysis for MHD systems studies. This report describes the effect 
of the progratmning details from a user point of view, but does not describe 
the Code in detail. 
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I „ INTRODUCTION 


The NASA-Lewis Chemical Equilbrium Code [1] had gained wide acceptance 
as a valuable and accurate tool for computing thermodynamic properties and 
chemical equilibrium compositions. In a companion publication [2], two of 
the authors have described a preprocessor program (GPSAP) written to 
facilitate general systems analysis studies. By adapting the NASA code to 
the format that GPSAP requires for component models, the very powerful 
systems analysis capabilities are obtained (e.g., execution time parameter 
alterations, automated sensitivity studies, and optimization features). 

Several additional capabilities were also added to the new code, to 
implement or simplify combustor studies, especially as related to MHD 
applications. These new features include: 

1) a routine for computing the electrical properties of the equili- 
brium mixture, 

2) a shorthand data input option for treating standard hydrocarbon 
fuel s , 

3) a facility for incorporating heat losses (as a fraction of thermal 
input) into the NASA HP mode computations, and 

4) a standard printed summary of equilibrium thermodynamic properties. 

Section II describes the modifications that have been made. Sec- 
tion III is a users guide, detailing the input data and giving information 
necessary for efficient systems studies. Two examples, both demonstrating 
the power of the modified code and illustrating its use are presented in 
Section IV. Throughout, it is assumed that the reader is familiar with 
both the original version of the NASA code and the fundamentals of GPSAP 


usage. 



II. CODE MODIFICATIONS 


A. GPSAP Compatibility 

As descr*ibed in Ref. 2, GPSAP requires certain standard features from 
the component models. Specifically, it is necessary that the data input, 
calculational, and data output modules be clearly separated entities, each 
with its own entry point, and that there are appropriate facilities for 
communicating data between the main program and the calculational block. 

Although it is quite generally possible to write PL/1 drivers properly 
modularized for FORTRAN models, the coupling is most easily accomplished 
with PL/1 models. The NASA code, as originally written, did not permit 
alteration of the input data at execution time. The most straightforward 
approach to remedying this situation was to translate into PL/1 the rou- 
tines for data input and reduction, output, and the calculational driver, 
and to rearrange the program logic according to GPSAP requirements. 
In doing so, the rocket, detonation, and constant volume equilibrium modes 
were discarded for the sake of execution efficiency. Aside from eliminating 
these options, the equilibrium calculations (EQLBRM, MATRIX, GAUSS, and 
CPHS) were laft unaltered in FORTRAN. 

The three PL/1 entry points to the new version are, 

NASAIN, 

NASA, 

NASAOUT. 

NASAIN simply reads the input data from file CARDIN and stores it. 
NASA performs any desired alterations to the input, reduces the resulting 
set of data to array form, calls EQLBRM repeatedly to perform the requested 
tasks, and, if instructed, prints various amounts of intermediate output. 
NASAOUT prints the summary of thermodynamic and electrical properties and, 
if instructed, stores the same on disk file FITDAT for post processing. 
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As In all GPSAP inodes, the different entry points are now completely 
independent. Thus, from a single call to MASAIN, NASA can be (and typically 
is) called repeatedly, the input data being altered each time from MAIN. 
The power of the approach derives from the fact that those input data 
alterations to one NASA call can (and typically do) derive from the results 
of a previous call. Thus, parameter optimizations may be performed in a 
single computer job. It has been extensively verified that, for a given set 
of input data, the answers obtained with the new code are identical to 
those obtained previously. 

B. Electrical Conductivity Procedure 
1 . Model Description 

The plasma electrical conductivity and mobility are computed by use of 

the algorithm of Demetriades and Argyropoulos [3]. This method calculates 

the electrical conductivity, using Grad's 13-moment approximation to derive 

a generalized Ohm's law for a multicomponent plasma. The electrical 

properties of the plasma are obtained in two successive approximations that 

are equivalent to the first and second approximations in the expansion of 

Sonine polynomials in the Chapman-Enskog approach. 

In this method, the coefficients of Ohm's law (electrical properties 

of plasma) are computed if three electron-neutral parameters are known as 

functions of electron temperature [4]: The effective collision cross 

( 2 ) 

section for momentum transfer, Qek* and the two weighting factors, A^ ^ 

and Ag^^, that characterize the electron-neutral interaction. 

The equations used in the calculation of the electrical properties of 
plasma that are required in Ohm's law are as follows: 

Ohm's law: E + U x B =(i]j + xJ x B 
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where the coefficients, a and x, are given by 


a - Oq/(1-a) 

X = (e ng)T^ 


where 
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or Agg = 8.7592-10® /n^^ 


These equations are the ones used In Ref. 3. However, In Ohm's law, 
the effective electric field term resulting from electron temperature and 
pressure gradients is neglected. Also, the effect of magnetic field on the 
Hall field coefficient (x) is neglected; it will be discussed later. 
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2. Input Data 

In the computation of the electrical properties of the plasma, the 
following species are considered: OH, Cs, Ar, CO, CO2, N2, O2, H2O, 
NO, KOH, SO2, K, and e (electrons). However, the NASA code considered 
all the possible species for the particular fuel used. The momentum cross 
sections, Q*s, and the weighting factors, A^2)'s and a( 5 )'s, for the 
specie considered are tabulated as functions of plasma temperature and are 
given in Table I. The values follow closely the data given by Spencer and 
Phelps. [5]. However, these data are subject to change should new experi- 
mental results become available. 

For the calculation of the electrical properties of plasma at other 
temperatures than those given in the table, linear interpolation is 
used. 
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Table 1 


0*s, k(2)*s and A(5)'s 

(as ftiiictions of temperature from 1700 to 3000K) 


OH 


17. E2, 
24. E2, 
87.00, 
61.91, 
-. 1981, 
1914, 
„921, 
.910, 


18. E2, 

25. E2, 

82.19, 

59.48, 

-.1969, 

-.1908, 

.920, 

. 917, 


1S.E2, 20. E2, 21. E2, 22. E2, 23. E2, 
26. E2, 27. E2, 28. E2, 29. E2, 30. E2, 
77.90, 74.05, 70.58, 67.42, 64.54, 
57.25, 55.18, 53.26, 51.48, 49.81, 
- .1950,- . 1 948,-. 1939,- . 1930,-. 1 922 
-.1991,-. 1895,-. 1 889, -. 1883,-. 1878 
.920, .920, .919, .919, .918, 

• t .917, .9^6, .916, .916, 


r 
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Cs 


17. E2, 
24. E2, 
685.48, 
458.56, 
-.2116, 
-. 3046, 
.501, 
.815, 


18. E2, 19. E2, 20. E2, 21. E2, 22. E2, 23 . E2 
25. E2, 26. E2, 27. E2, 28. E2, 29. E2, 30. B2 
645.31, 607.81, 572.99, 540.78, 511.07, 

435.45, 414.24, 394.77, 376.91, 360.50, 

-.2334, -.2520, -.2677, -.2806, -.2909, 

-.308, -.3100, -.3101, -.3085, -.3055, 

.539, .582, .626, .673, .720, 

.862, .907, .951, .994, 1.034, 


483.71, 

345.44, 

-.2989, 

-.3012, 

.768, 

1.073, 


17. E2, 

18. E2, 

19.E2, 

24. E2, 

25. E2, 

26. E2, 

.494, 

.526, 

.560, 

.745p 

.783, 

. 821, 

.6356, 

.6574, 

.6723, 

.6919, 

. 6904, 

.6880, 

3.085, 

3 .062, 

3.033, 

2.863, 

2.832, 

2.807, 


20. E2, 

21. E2, 

22. E2 

27. E 2, 

28. E2, 

29. E2 

.596, 

. 632, 

. 669, 

.860, 

.899, 

. 937, 

. 6822, 

.6882, 

. 6913 

. 6851, 

. 6819, 

. 6784 

3. COO, 

2.965, 

2.930 

2.775, 

2.750, 

2.727 


23. E2, 

30. E2, 

.707, 

.976, 

.6924, 

.6749, 

2.896, 

2.705, 


CO 


17. E2, 
24. E2, 
13.050, 
15.219, 
.3617, 
.3933, 
1.771 , 
1.829, 


18. E2, 19. E2, 20. E2, 21. 
25. E2, 26. E2, 27. E2, 28. 
13.360, 13.669, 13.979, 
15.523, 15.823, 16.117, 
.3668, .3722, . 3777, 

.3948, .3951, .3943, 

1.800, 1.230, 1.840, 

1.807, 1.779, 1.746, 


E2, 22. E2, 23. E2, 

E2, 29. E2, 30. E2, 

14.290, 14.60 1, 14.911, 
16.403, 16.680, 16.947, 
.3823, .3872, .3908, 

.3924, .3894, .3855, 

1.850, 1.851 , 1.844, 

1.-08, 1.668, 1.626, 


CO 2 

17. E2, 
24. E2, 
19. 15, 
13.77, 
-. 1802, 
-.1796, 
.857, 

. 957, 


18. E2, 19. E2, 20 
25. E2, 26. E?, 27 
18.13, 17.22, 

13.25, 12.77, 

-.1823, -.1836, 
-.1769, -.1738, 
.870, .883, 

.974, .992, 


E2, 21, E2, 22. 
E2, 28. E2, 29. 
16. 39, 15.64, 

12.33, 11.92, 

-.1841, -.1840, 
-.1700, -.1658, 
.897, .911, 

1.010, 1.029, 


E2, 23. E2, 

E2, 30. E2, 

14.96, 14.34, 

11.55, 11.20, 

-.1831, -.1816, 
-.1610, -.1558, 
.926, .942, 

1.049, 1.070, 
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Table 1, Cont'd 


u 


17. E2, 

18. B2, 

19. B2, 

20. E2, 

21. E2, 

22. E2, 

23. E2, 

24. B2, 

25. E2, 

26. E2, 

27. E2, 

28. E2, 

29. F2, 

30. B2, 

9.980, 

9.178, 

9.279, 

9.358, 

9.445, 

9.532, 

9.620, 

9.710, 

9.804, 

9.903, 

10.007 

,10. 115, 

10.229 

,10.349 

.2757, 

.2741, 

.2737, 

.2747, 

.2769, 

. 2004, 

.2852, 

.2911, 

.2980, 

.3058, 

.3143, 

.3234, 

. 3330, 

.3427, 

1.428, 

1.444, 

1.467, 

1.495, 

1.529, 

1 .568, 

1.611, 

1 657^ 

1 .705, 

1.753, 

1.801, 

1.847, 

1.891, 

1.931, 

17.E2r 

18. E2, 

19. E2, 

20. E2, 

21. E2, 

22. E2, 

23. E2, 

24. B2, 

25. E2, 

26. E2, 

27. E2, 

28. E2, 

29. E2, 

30. S2, 

n.624. 

5.696, 

5.765, 

5.832, 

5.896, 

5.958, 

6.018, 

6.075, 

6.129, 

6.181, 

6.231, 

6.278, 

6.322, 

6.364, 

.2896, 

.2896, 

.2697, 

.2898, 

.2898, 

.2895, 

.2890, 

. 2881, 

. 2869, 

.2855, 

.2837, 

.2816, 

.2793, 

.2768, 

1.498, 

1.501, 

1.502, 

1.501, 

1.496, 

1.490, 

1.482, 

1.472, 

1 .460, 

1.447, 

1.433, 

1.419, 

1.404, 

1.388, 

17. E2, 

18. E2, 

19. E2, 

20. E2, 

21. E2, 

22. £2, 

23 ,E2, 

24. E2, 

25. E2, 

26. E2, 

27. E2, 

28. E2, 

29.E2. 

30. E2, 


Ti08. 53?,, 102. 550,97.229, 92.443, 88.120, 84.197,' 80.620, 
77.344, 74.333, 71.555, 68.985, 66.598, 64.377, 62.305, 

-. 1968, -.1955, -. 19tv2, -. 1931, -.1921, -.1911, -.1903, 

-.1895, -.1887, -.1881, -.1874, -.1869, -.1863, -.1859, 

.923, .923, .922, .922, .921, .920, .919, 

.919, .918, .917, .916, .916, .915, .915, 
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SjQ2 

17. E2, 18. E2, 19. E2, 20. E2, 21. P2, 22. F2, 23. E2, 

24. E2, 25. F2, 26. E2, 27. E2, 28. E2, 29. E2, 30. E2, 

84. 93, 80.25, 76.08, 72.33, 68.95, 65.88, 63.08, 

60.51, 58.15, 55.98, 53.97, 52.098, 50.359, 48.736, 

-.197, -.1957, -.1945, -.1934, -.1923, -.1914, -.1905, 

-.1897, -.1890, -.1884, -.1878, -.1872, -.1867, -.1863, 

.923, .923, .922, .922, .921, .920, .919, 

.918, .917, .916, .916, .915, .915, .914, 

K 

17. E2, 18. E2, 19. E2, 20. E2, 21. E2, 22. E2, 23. E2, 

24. E2, 25. E2, 26. E2, 27. E2, 28. E2, 29. E2, 30. E2, 

300.143, 284.791, 271.027, 258.609, 247.345, 237.076, 227.674, 
219.029, 211.053, 203.669, 196.814, 190.432, 184.477, 1', IQB , 
-.1679, -.1669, -.1661, -.1654, -.1648, -.1644, -.1«40, 

-.1636, -.1633, -.1630, -.1627, -.1624, -.1620, -.1616, 

.921, .919, .917, .915, .914, .912, .911, 

.911, .911, .911, .912, .913, .915, .917, 



3. Calculation of the Hall Parameter, 3 

According to the algorithm adapted [4] for the calculation of electrical 
properties of plasma, the Hall parameter, 3, is 
3 = X-a-B, 

where B is the magnetic field. 
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However, x is also a function of the magnetic field, where [4] 


X = (en^) 


and 






)■ 


eB 


e m 


The above equation for x can be simplified, after substituting for the 
* 

variables A, t , and oig, and the equation becomes: 


xMenJ-' 

32 


where 


^0 " 

Cl . y (-|) v/. 

1 0 I ^ e 0 


The influence of the magnetic field on the Hall coefficient, x, has 
been investigated for zero magnetic field, up to 20 T for a typical MHD 
plasma composition. The results are presented in Fig. II-l. These results 
indicate that there is a very weak dependence of the coefficient, X, on the 
magnetic field. This is a result of and C.| being very small and, also* 
almost equal . 

Based on these results, the effect of the magnetic field on the Hall 
coefficient, x, is neglected, and the equation for the computation of X is 
reduced to 


X = (eng)T^ 

4. Error Bounds on Electrical Conductivity 

For the same method used to compute the electrical conductivity of the 
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Halt coefflc Lents at Cohen points 


Fig. II-l 
Effect of the 
magnetic field 
on the Hall 
Coefficient (as 
defined in 


Ref. 4) 


a 



plasma, there are two major sources of errors that can contribute to the 
uncertainty of the computed values for conductivity. The first source of 
error is the thermochemical data for the heat of formation of some of the 
species, which can significantly affect the concentration of the electrons 
in the plasma. The most important species are KOH, A102» and CO 2 . The sec- 
ond source of error is the values for the momentum collision cross sections. 
Reference 7 discussed the error bounds on the electrical conductivity, which 
result from the uncertainty of the data of the first kind. At ANL, indepen- 
dent calculations were done to Investigate such effects. 
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For CO 2 : Using the JANAF data of December, 1975, a is reduced by 5 * at 

a reference point (p = 3.5 atm, T = 2600 K, for Montana Rosebud 
coal with 5 % moisture, 85 % slag rejection, and 34 % oxygen 
enrichment) as compared to the case where CO 2 is neglected. 
Based on the investigation of References 7 through 9, it has been 
shown that, with the new data for the heat of formation for 
CO 2 , there is very little effect of CO 2 on the electron con- 
centration in the plasma and, correspondingly, very little 

}>. 

effect on the electrical conductivity of the plasma. On this 
basis, until the JANAF tables are officially updated, CO 2 is 
neglected in our analysis. 

For AIO 2 : Using the old JANAF data of December, 1968, a is reduced by 

approximately 13 %. The error bounds on the heat of formation 

in the new data (12/75) are smaller, which makes the range of 

uncertainty smaller ( 'v-13 %) than before. 

For KOH: Reference 7 showed that s can vary by as much as 50 % because 

of the uncertainty of 3 kcal/mol in the heat of formation. 

Therefore, it is important to reduce experimentally the range 

of uncertainty. In the meantime, ANL uses the JANAF medium 

value. 

C. Heat Losses 

When run in the HP mode, the original NASA code yielded a single 
calculational point for each pressure input. The heat of formation of each 
reactant, as input or as computed, was multipled by the mole fraction of 
that reactant, and the products were summed to give the total enthalpy, 
HSUBO. Chemical equilibrium was then calculated for that total enthalpy 
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and each specified pressure. For many purposes, especially combustor 
studies, it is valuable to be able to find equilibrium properties as a 
function of combustion chamber heat loss. The new version has been ad- 
justed to allow these studies. 

These heat loss studies are performed by computing chemical equili- 
brium at an enthalpy that is reduced by some amount from HSUBO. Because 
enthalpy is a relative measure, a reference value, REFH, was added to the 
inputs. This is typically the enthalpy of the same reactants when brought 
to equilibrium at some standard set of conditions, and is treated as the 
enthalpy when heat loss is '100 For a given fractional heat loss, f, 
the enthalpy at which equilibrum is computed is 

H=HSUB0 - f (HSUBO - REFH). 

To save space, the heat loss fractions are input and stored in the 
T array, the same array as is used for temperatures in TP mode calculations. 
Because this array is internally counted by looking for the first zero 
element, and because no heat loss should be an option, the values actually 
input to T are from (1-f), rather than f. 

Thus, for example, to compute the combustor temperature for a 2.5% 
combustor heat loss (based on the total thermal input), the HP input mode 
is used and the value 0.975 is entered in the T array of the NASA input. 

D. Shorthand Data Input 

The reactants for virtually all MHD combustor studies, and for a large 
number of other applications, conform to a standard model. A hydrocarbon 
fuel (coal, oil, natural gas, or some combination of them) of given chemical 
analysis and heating value is burned in a preheated air/oxygen mixture with 
some weight fraction of potassium, injected as K 2 CO 3 . The air, fuel, and 
K 2 CO 3 may be accompanied by some given amount of water. 
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Although, relatively straightforward, inputting this standard 
analysis data by REACTANTS cards is rather tedious. For use in coal 
combustion, the weight fractions must be adjusted from the 'as received' 
to the 'as determined' basis, in order to correct for various moisture 
content and ash rejection ratios. The weig'it fraction of air and pure 
oxygen must be adjusted to meet the desired 02-enrichment and stoichio- 
metric ratios. The weight fraction of K2CO3 must be adjusted to the 
proper fraction of potassium seeding. And a fictitious heat of formation 
for one reactant is required, so that enthalpy of the total mixture will 
correspond to the given heating value of the fuel as correctly adjusted to 
the 'as determined' basis. 

To simplify treatment of this standard problem, a new input group, 
FUEL, has been added. As described in detail in Section IV, a FUEL 
group defines the following: 

the chemical composition of the fuel, 
fuel heating value, 
moisture content, 
ash rejection, 

stoichiometry of the combustion, 

02-enrichment, 
fuel temperature, 
air temperature, 

O2 temperature, 

potassium (seed) weight fraction, 
water-to-seed weight fraction. 

Within the computational entry, NASA, the relative weights and the 
fictitious enthalpy of the carbon are computed. All the information from 
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the FUEL group is then stored in the appropriate internal arrays - duplica- 
ting the situation that would have occurred if the reduced data had been 
read from REACTANTS cards. 

The computations that are performed in this setup process are as 
follows: 

1) From the ultimate analysis of the fuel, the total fuel weight 
(T0T_WT) is readjusted to yield the specified value of fuel moisture 
content (FUEL .H20_ADJUST) . 

T0T_WT = T0T_WT*(1 .0-H20_FRAC_FUEL)/(l .0-FUEL. H20_ADJUST) 

The 'as received' heating value for the fuel (FUEL.HHV) is readjusted to 
account for this moisture change. 

HHV = HHV*(1 .0-FUEL.H20_ADJUST)/(l .0-H20_FRACJ"UEL) 

2) From the specified fraction of ash rejection (FUEL.ASH_REJECT) , 
the weights, weight fractions, and HHV are readjusted. 

SLAG_WT = total weight of ash content as received, 

HHV = HHV*TOT_WT/(TOT_WT-ASH_REJECT*SLAG_WT), 

T0T_WT = TOT_WT-ASH_REJECT*SLAG_WT, 

X = X*{1 .0-ASH_REJECT), 

where X is the weight of any ash species (e.g., Si02, SO3). 

3) From the input carbon (C), hydrogen (H), sulfur (S), and oxygen 
(0) content of the fuel and the stoichimetry (STOICH), the total weight of 
oxygen (OWT) needed in the oxidizer mixture is determined. 

OWT = ST0ICH*(32.0*(C/12.0111 ) + 0.5*(H/2.016) + (S/32.066) )-0, 
where the molecular weights of O2, C, H2, and S are 32.0, 12.011, 

2.016, and 32.066, in that order. 

4) From the oxygen enrichment, the total amount of air (AIRWT) is 
determi ned . 

AIRT = OWT/02 AIR 
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5) From the air weight and oxygen weight, the weight of pure oxygen 
(02WT) is determined, and, from the humidity (H20A), the amount of water 
(H20WTA) in the oxidizer mixture is determined. 

02WT = 0WT-0.2315*AIRWT 
H20WTA = H20A*AIRWT, 

where the mass ratio of oxygen in dry air is 0.2315. 

6) From the air temperature (AIRT), its enthalpy (AIRH) is obtained 
from the polynomial fit. 

AIRH = -1.174E-10*(AIRT)**4+3.823E-7*(AIRT)**3 
+2.3485E-4*(AIRT)**2+6.713*(AIRT)-2059.76 

7) From the seed fraction (SF), the fuel weight (FW), and the mois- 
ture fraction (H20S), the amount of K2CO3 (SEEDWT) and the amount of 
water (H20WTS) accompanying it are determined. 

SEEDWT = SF*(FW+02WT+AIRWT+H20WTA)/(0.5658-SF*(l+H20S)) 

H20WTS = H20S*SEEDWT, 

where 0.5658 is the mass ratio of K in K2CO3. 

8) From the heating value (HHV) of the fuel, its total weight, and 
its carbon, hydrogen, and sulfur content, the appropriate fictitious carbon 
heat of formation ('Hf') is determined. 

Hf = -94052+12. 011/C*(HHV*FW-H/2. 016*6837. 0-5/32.066*70947.0), 
where the heat of formation values for CO2, H2O, and SO2 are -94052, 
68317, and 70947 cal/g, in that order. 

These steps are applied to each FUEL group automatically, loading the 
internal arrays within the NASA code as if a set of REACTANTS cards had been 
furnished. 
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E. Thermodynamic and Electrical Properties Summary Output 

When the output routine, NASAOUT, is called, it causes the printing 
of a summary sheet containing the overall thermodynamic properties of the 
equilibrium mixture. This summary consists of the following labeled columns 
of data. 


Col umn Label 


Description 


1 

2 

3 

4 

5 
5 

7 

8 
9 


P 

T 

MW 

H 

S 

RHO 

CP 

SON VEL 
SIG 


11 


NE 


pressure in atmospheres (input quantity) 
Temperature in Kelvins (input quantity in TP mode) 
Equilibrium molecular weight 
Enthalpy in joules/kg 

Entropy in joules/kg *K (input quantity in SP mode) 

3 

Density in kg/m 

Specific heat in joules/kg’K 

Sonic velocity in m/s 

Electrical conductivity in S/m (present only if 

conductivity calculation required) 

20 3 

Electron number density in 10 /m (present 

only if conductivity calculation is required) 


These numbers are stored internally as the rows of the PL/1 -based struc- 
ture array, NASAP->OUT.A (cf Ref. 2) whenever the calculational entry, NASA, 
is called. Thus, for example, the values 


NASAP->0UT.A(10,J) 

for J = 1, 2, . . . , NASAP->IFIN refer, in order, to the conductivities 
resulting from each equilibrium point calculated. Because the pointer, 
NASAP, is an external variable, all of these data points are known to (and 
may be used in) the MAIN Program. 


19 



III. USER'S GUIDE 


This section will describe in detail the input data to the NASA code 
and the output that may be obtained with the various options. The variable 
names into which the input data is stored are also listed and explained, 
because knowledge of these names is necessary to effect control from GPSAP 
over parameters in NASA. 

A. Input Data 

Data input to NASAIN, contained in file CARDIN, is divided into cate- 
gories, each introduced by a keyword. The keywords and the type of input 
data following each a»'e: 



Keyword 

Type of data 

1) 

rEACTANTS 

Description of chemical reactants introduced into the 
combustor 

2) 

FUEL 

Shorthand description of reactants 

3) 

NAMELIST 

Type of problem, prescription of "givens," selection of 
options 

4) 

OMIT 

Chemical species to be excluded in computing equilibria 

5) 

INSERT 

Condensed species to be included in first iteration 

6) 

END 

Signifies the completion of problem prescription 


These keywords must appear, starting in Column 1. Any additional data 
on the same line (card) with REACTANTS, FUEL, NAMELIST, or END key cards is 
ignored; thus, comments may be entered in those lines. 

REACTANTS, OMIT, and INSERT all trigger 'GET EDIT' statements, and the 
data must, therefore, occur in specified columns. 

FUEL and NAMELIST trigger 'GET DATA' statements, allowing more free- 
form input. 
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END returns program control to the calling (MAIN) program. Inclusion 
of this key card In the program vocabulary allows the other Input categories 
to occur In any order. The various categories are discussed In detail below. 
1 . REACTANTS 


Except for omitting the division of reactants Into "fuels" and "oxi- 
dants," the REACTANTS group format has not been changed from that of the 
original NASA-Lewls code. Thus, the REACTANTS groups from any existing 
data sets for that code can be used directly In the present one. The 
names, meanings, and Input format for each REACTANTS card (from 1 to NREAC, 
where NREAC Is the total number of reactants read) are as follows: 

Variable Meaning Format Col umns 

NAME(N,I),NUM(N,I) Atomic symbol and formula 5(A(2) ,F(7,5) ) 1-45 

numbers for element I 
(maximum 5) of reactant N 
(If NAME(N,5) ='00,' the 
enthalpy Is calculated by 
the code). 


PECWTS(N) Relative weight or number F(7,5) 46-52 

of moles of reactant N In 
total reactants. Arbitrary 
base. 


MOLE If M0LE='M' for any N, then A(l) 53 

variable MOLES Is set to 
'I'b, signifying that aJl 
PECWTS are Interpreted as 
relative numbers of moles. 

MOLES Is alterable at 
execution time, MOLE Is not. 


ENTH(N) Assigned enthalpy In F(9,5) 54-62 

cal/g*mol [see note 
regarding NAME(N,5)]. 


FAZ(N) 


'S', If solid A(l) 63 

■L', If liquid 
'6', If gas 
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RTEMP(N) 

Reactant temperature in K 

F(8,5) 

64-71 

DENS(N) 

(optional) reactant density 

F(8,5) 

73-80 


in g/cm^ 




Note that column 72 is ignored. All of the above, except MOLE (but 
including MOLES), can be altered with 'AAA' statements at execution time, as 
shown in the examples in Section IV. 

2. FUEL 

The FUEL keyword is a local extension, designed to simplify data input 
for the most common MHO applications. After the keyword, any members of 
the structure FUEL are input as free-form assignment statements (e.g., 
FUEL.O = 86.4) separated by commas and/or any number of blanks. The end 
of a FUEL group is designated by a semicolon. The names, meanings, and 
default values of the variables input in a FUEL group are 


Vari 

able 

Meaning 





Def au 

FUEL 

.C 

Relative 

weight 

of 

carbon 

in fuel 

0.0 

FUEL 

.H 

Rel ative 

weight 

of 

hydrogen in fuel 

0.0 

FUEL 

.S 

Relative 

weight 

of 

sulfur 

in fuel 

0.0 

FUEL 

.N 

Relative 

weight 

of 

nitrogen in fuel 

0.0 

FUEL 

.0 

Relative 

weight 

of 

oxygen 

in fuel 

0.0 

FUEL 

.H20 

Relative 

weight 

of 

H^O in 

fuel 

0.0 

FUEL 

.Si02 

Relative 

weight 

of 

SiO^ ii 

T fuel 

0.0 

FUEL 

.S03 

Relative 

weight 

of 

SO^ in 

fuel 

0.0 

FUEL 

.AL203 

Relative 

weight 

of 

AI 2 O 3 

in fuel 

0.0 

FUEL 

.CAO 

Relative 

weight 

of 

CaO in 

fuel 

0.0 
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FUEL.MGO 

Relative weight of MgO in fuel 

0.0 

FUEL.FE203 

Relative weight of Fe202 in fuel 

0.0 

FUEL.K20 

Relative weight of K2O in fuel 

0.0 

FUEL.TOT_WT 

Total weight to which the above 
components are normalized 

100.0 

FUEL.T 

Temperature of fuel 

300.0 

FUEL.HHV 

Higher heating value of fuel in cal/g 

0.0 

FUEL .DENS 

density of fuel 

0.0 

FUEL.STOICH 

Ratio of total oxygen supplied in 
fuel, air, and pure 0^ to total 
oxygen required for complete 
combustion of all carbon, hydrogen, 
and sulfur 

1.0 

FUEL.02_AIR 

Weight ratio of the total O2 to the 
air + O2 mixture supplied^ 

0.2315 

FUEL-H20_AIR 

Weight of water carried by the air, 
as a fraction of the dry air weight 

0.0 

FUEL.AIR_TEMP 

Temperature of air and its water 

300.0 

FUEL.02_TEMP 

Temperature of the pure O2 

300.0 

FUEL.SEED_FRAC 

Weight fraction of potassium in 
total to be introduced as K2CO2 

0.01 

FUEL.H20_SEED 

Weight of water accompanying the 
K2CO3 as a fraction of the 
K2CO3 weight 

0.0 

FUEL.H20_ADJUST 

Weight fraction of H2O required in 
the fuel. Thus, if FUEL.H20 ADJUST 
is not equal to the moisture frac- 
tion in the fuel, the weight fractions 
of water and the HHV of the fuel are 
readjusted to this new moisture frac- 
ti on . 

0.0 

FUEL.ASH_REJECT 

Fraction of the total ash content of 
the fuel to be rejected 

0.0 


Notes: 

(1) If FUEL.C is not v^ntered, entire FUEL group is ignored. 
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(2) The seed fraction is calculated as a fraction of the total of the 
reactants generated by this FUEL group. 

(3) As many as three FUEL groups and one REACTANTS group may be 
entered to describe the reactants for a given problem. At input, 
each FUEL group is tagged by the keyword, FUEL, and input as 

FUEL.xxx=vvv 

(4) For modification at execution time, members of the second and 
third FUEL groups entered are modified by 

AAA 'FUEL2.XXX' vvv 

and 

AAA 'FUELS. XXX' vvv, 
respectively. 

3. NAMELIST 


The keyword, NAMELIST, triggers a GET DATA statement, which reads a 
number of parameters and switches that control the calculations to be 
performed. Although the basic idea is the same as in the original NASA- 
Lewis code, considerable modification of this input categroy was necessary, 
both to accomodate the PL/I environment and to serve the locally available 


options. 

Variable 

The variables, types, meanini 
Dimension Type 

TP 

1 

BIT(l) 

If 

•I'B, 




at 

given 

HP 

1 

BIT(l) 

If 

'I'B, 




at 

given 

SP 

1' 

BIT(l) 

If 

'I'B, 




at 

given 

OTTO 

1 

BIT(l) 

If 

'O'B, 


Meaning 


Default 


'O'B 


JANAF chemical table 
search is skippod whenever the 
preceding NAl^A call had the same 
elements and omissions 


'O'B 

'O'B 

'O'B 
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Variable Dimension Type Meaning Default 




lU 



MMHG 

1 

BIT(l) 

If 'I'B, pressures are in mm Hg 

'O'B 

INGH 

1 

BIT(l) 

If 'I'B, pressures are in in. Hg 

'O'B 

PS I A 

1 

BIT(l) 

If 'I'B, pressures are in psi 

'O'B 

NSQM 

1 

BIT(l) 

If 'I'B, pressures are in 
N/m2. Note that, if none 
of MMHG, INHG, PSIA, or NSQM, 
then pressures are in atmospheres 

'O'B 

FIT 

1 

BIT(l) 

If 'I'B, summary output is 
written onto file FITDAT for 
subsequent polynomial fit 

'O'B 

CONDUC 

1 

BIT(l) 

If 'I'B, conductivity routine 
is called after equilibrium 
calcul ation 

'O'B 

IONS 

1 

BIT(l) 

If 'I'B, ionized species are 
included in equilibrium 
calculation 

'I'B 

PRINT 

1 

BIN(15) 

Determines the amount of inter- 
mediate output 

0 

IDEBUG 

1 

BVIN(31) 

If > 0, causes equilibrium 
iterations to be output 

0 

TRACE 

1 

DEC(6) 

If > 0, signifies that, when 
mole fractions are printed, all 
species with fraction > TRACE 
are shown 

0 

REFH 

1 

DEC(6) 

Reference enthalpy: enthalpy 

of the same reactants if brought 
to equilibrium at standard temp 
and pres. In TP mode, this is 
used as a "zero set" for the 
reported enthalpies. In HP mode, 
it is used as a reference point 
to define heat loss 

0 

SO 

1 

DEC(6) 

Total entropy at which equilibria 
are calculated in SP mode 

0 

P 

26 

DEC(6) 

Schedule of pressures at which 

26*0 


equilibria are to be calculated. 
End of the schedule is signified 
by the first 0 
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Variable Dimension Type 


Meaning 


Default 


T 26 DEC(6) If TP mode, schedule of temperatures 26*0 

(in K) at which equilibria are 
calculated. If HP mode, schedule 
of heat retention fractions at which 
equilibria are calculated (see below) 

The use of the T array in HP mode is a locally written extension to 
the original NASA-Lewis code. In this application, the values 1 - T(I) 
represent the combustor heat loss as a fraction of the enthalpy difference 
between the reactants as introduced to the combustor and the enthalpy 
REFH. Thus, e.g., if REFH is the enthalpy of the same reactants at 300 K, 1 
atm, then a value T = .96 represents a heat loss of 4% of the total 
thermal input. 

Note that, in a NAMELIST input from a call to NASAIN, all variables 
are redefined to, their default values before the NAMELIST is read. This 
can be important when GPSAP is used to do parameter sweeps and optimizations. 
4. OMIT and INSERT 

OMIT and INSERT cards are identical in both function and format to 
those in the original NASA-Lewis code. The keyword is followed by the 
names (chemical symbols) of from one to four species, starting in columns 
16, 31, 46, and 61. The names on OMIT cards signify particular species that 
are to be deleted from consideration in computing the equilibrium composi- 
tion. INSERT cards contain the names of condensed species that are assumed 
to be present when the program is begun. The variables, NSERT and NOMIT, 
contain, respectively, the number of inserted and omitted species; ENSERT 
(3,30) and 0MIT(3,3) [both CHAR(4)] contain the species names. At execution 
time, the last species in either list may be removed by decreasing NSERT or 
NOMIT. A species may be changed to a new one by assigning the appropriate 
new character strings to the array elements. The name of the Nth species 

is left- justified in the 12-character field, 
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(XXX{1,N). XXX(2,N), XXX(3,N)), 

where XXX stands for either OMIT or ENSERT. If the old species name is 
longer than the new one, the extra spaces must be explicitly blanked. A 
new species may be added to the list by increasing NSERT or NOMIT and 
assigning the new name to the array. 

5.0 END 

The END card does not signal additional inout but, rather, indicates 
that a problem is now fully specified and that control should be returned 
to the main program for further processing. 

B. Program Output 

Primary output from the calculation is the thermodynamic and electrical 
properties summary printed upon a call to NASAOUT. This output consists of 
a table for all the points on the P and T schedules of the following: 
pressures (atm) 
temperatures (K) 
molecular weights (kg/kg-mol) 
enthalpies(J/kg) 
entropies (J/kg*K) 

3 

densities (kg/m ) 
specific heats (J/kg-K) 
gas constants 
sonic velocities (m/s) 
condUiCtivities (S/m) 

20 3 

electron number densities (10 /m ) 

These same values are also accessible from the driver as NASP”>OUT(I,J), 
where I goes from 1 to 13 and refers to the properties (in the order given) 


and J goes from 1 to the number of output points (T varies more rapidly 
than P). 

In addition, various amounts of intermediate output may be obtained, 
depending on the value of PRINT. 

If PRINT =0, all that is obtained is a copy of the input data and any 
error messages. Error messages will be on file FT06F001 if errors occurred 
in the equilibrium calculation and on SYSPRINT otherwise. The messages and 
their meanings are identical to those in the original code and will not be 
repeated here. 

If PRINT=1, the output also includes a list of the species considered 
each time the JANAF thermodynamic data are scanned, the value (HSUBO) of 
the total enthalpy of the reactants as introduced, the number of gram atoms 
of each element present per kg of the mixture (BO's), a listing of any 

inserted species, and a list of the REACTANTS cards that are equivalent to 

any FUEL groups used and a list of the reactants, weight fractions (normal- 
ized to a total of 1), enthalpies (as input or calculated), phases, refer- 
ence temperatures, and densities (if input), as well as the net mixture 
density at input. 

PRINT=2 adds a brief summary of the net thermal and electrical proper- 
ties printed once every 13 points on the calculation schedule. 

PRINT =3 yields all of the above plus 1) a listing of the mole fractions 
for every species of mole fraction greater than or equal to 5 X 10~® if 
TRACE=0 (TRACE otherwise), 2) a trace of the LaGrange multiplier values, 
and 3) the names of added or removed condensed species as the equilibrium 
calculation proceeds (this last is on file FT06F001). 
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IV. ILLUSTRATIVE EXAMPLE 


Two examples will both demonstrate the method of using the new GPSAP- 
driven NASA code and illustrate the new capabilities it affords. Any 
confusion as to variable names and meanings may be resolved by reference to 
Section III. 


A. Example 1 


1. Problem 

Determine the thermodynamic and electrical propeties of a coal burned 
in an air/oxygen mixture. The coal has an 'as received' composition 
analysis of: 

Species Weight 


C 

H 

S 

N 

0 

Si02 

SO3 

A1 2O3 

CaO 

MgO 

Fe203 

H2O 


52.2 
3.38 
0.8 
0.8 

11.02 
1.732 
2.055 
1.214 
1.481 
0.517 
0.499 

24.30 


Before being burned, the coal is dryed to 5% moisture content; further- 
more, it is assumed that the ash rejection is 75%. (For the analysis, it 
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iis assumed sufficient to remove 75% of the ash species before the combus- 
tion, rather than 75% of the condensed species during the combustion process.) 

The oxygen/air mixture is to have a preheat temperature of 1000 K and 
30 mass % O2. Sufficient air will be added to obtain a stoichiometric 
ratio of 0.92. The air is assumed to have moisture content of 1 mass %. 

The combustion process is also seeded with 1 mass % potassium. The 
seed and coal are input at a temperature of 300 K. Also, the higher 
heating value of the coal is assumed to be 4968 cal/g. 

Specifically, the flame temperatures of the combustion are to be deter- 
mined for pressures is of 4, 5, 6, 7, and 8 atm with no heat loss. In addi- 
tion, the thermodynamic and electrical properties are to be determined for 
pressures in the range of 1-8 atm, in 1 atm increments, and for temperatures 
in the range of 1600-3000 K, in 200 K increments. 

2. Input 

The total input required for this analysis is shown in Fig. IV-1 . 
This input will now be discussed in detail, on the basis of the line 
numbers shown at the left in the figure. 

Lines 1 to 29 comprise the standard JCL, common to all GPSAP analysis. 

The in-line procedure, SYSTEM, that is executed (line 29) is composed of three 
steps: ONE, PLO, and TWO. Step ONE executes GPSAP, which reads the 

structure in file STRUCIN and expands it into a PL/1 MAIN program on scratch 
file SYSDRV. Step PLO compiles that program, and step TWO loads it, using 
library SYSLTB. 

Lines 33 to 38 complete the JCL, defining the various data files 
referenced by NASA. CNDCT contains the table of coefficients used to 
compute the conductivity. JANAF contains the thermochemical data (JANAF 
table) used by NASA. FITDAT is a scratch file used to store the NASA 
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Figure lv-1 


I 


1 . JOB 

2 . CDA 

3. //SISTEH PPOC 

4. //OHE EXEC PGH=GPSXP 

5. //STEPLIB DD DS»=C 11 5. B25236 . S? SLIB . LORD ,DISP=SHB 

6. //STBDCT DD DDHXHE=STBOCIN 

7. //SYSPBINT DD S YSO DT=X , DCB= (BE2 F M=VB A, LBECL= 1 37 ,BLKS IZ E= 1 51 1) 

8. //SYSDBV DD DISP= (NEW, PASS) ,DMIT=S ASCB, SPACE= (TBK, (2, 1) ) 

9. //PLO EXEC PGH=IELOAA,PABB='NS,NA,NX,HAG,HOESD,HSTG* 

10. //STEPLIB DD DSK=PLI.0PT.L1NKLIB,DISP=SHB 

11. //SYSIH DD DSN=*. ONE. SYSDBV, DISP= (OID, DELETE) 

12. //SYSLIH DD ONTT=SASCB,SPACE= (CYL,4) , 

13. // DISP = (HEW,PASS) ,DCB= (PECFM-PB,LBECL=80,BLKSIZE=3120) 


14. 

//SYSPBINT 

DD 

DDHMY 

15. 

//SYSPDHCH 

DD 

DDMHY 

16. 

//SYS DTI 

DD 

SPACE= (CYL,4) ,DNIT=(SASCE) 


17. //TWO EXEC PGM=LOADEB,BEGION=150K,COND=(9,LT,PLO) 


18. 

//SYSLIB 

DD 

DSN=SYS1.PLIBASE,DISP=SHB 

19. 

// 

DD 

DSN=SYS 1 . FORT LIB, DISP=SHP 

20. 

// 

DD 

DSN=C115.B25236.L0AD,DISP=SHB 

21 . 

//SYSLIN 

DD 

DSN=*. PLO. SYSLIN, DISP= (OLD, DELETE) 

22. 

//SYSLODT 

DD 

SYSODT=A,DCB= (EECFH=FB,LPECL=121 , BLKSIZE=1573) 

23. 

//SYSPNCH 

DD 

DOMMY 

24. 

//SYSPBINT 

DD 

SYEOOT=A,DCB= (BECFH=VBA, IBECL = 1 37, BLK SIZE= 1 51 1) 

25. 

//SYSPDNCH 

DD 

DOHMY 

26. 

//PBOGBEP 

DD 

SYSODT=A,DCB= (RECFH=FB, LBPCL = 133,BLKSIZE=1596) 

27. 

//FT06F001 

DD 

SYS00T=A,DCB= (BECFS=FB,LBECL=133,BLKSIZE=1596) 

28. 

// “END 



29. 

// EXEC SYSTFB 

30. 

//ONE.STBDCIN 

DD ♦ 

31. 

HNN NASAIN 

NASA NASAOUT ; 

32. 

NNN NASAIN 

NASA NASAODT ; 


33. //TWO.CNDCT DD DSN=B25445. TABLE . DATA, DISP=SHB 

34. //JANAF DD DSN=B25445.THEBMO.DATA,DISP=SHB 

35. //SUHHAHY DD SYSOUT=A,DCB= (BECFM=FB,LBBCL=133,BLKSIZE=1596) 

36 //FITDAT DD DISP= (NEW,PASS) , aNIT=SASCB, SPACE= (TBK, (5, 5) ) , 

37 // DCB=(BECFH=FB,LBECL=80,BLKSIZE=1520) 

38 //CABDIH DD * 

39 FUEL 

40 F0BL.C=52.2 FDFL.H=3.38 FDEL.S=0.8 FUEL. 11=0. 8 POEL.O=11.02 

41 FDEL.SI02=1 .732 FOEL. S03=2. 055 FDEL. AL203=1. 214 F0EL.K2O=0.0 

42 FOEL. CAO=1. 481 FDEL .HGO=0 . 51 7 FDEL .FE203=0 . 499 FDEL. H20=24. 3 

43 FDEL. ASH BEJECT=0.75 FDEL.H20 ADJDST=0.05 

44 FDEL. SEED FBAC=0.01 FDEL.H20 SEED=0.00 

45 FDEL.HHV=4968. FDEL.STOICH=0792 

46 FDEL.T=300. FDEL.AIB TENP=1000. FDEL. 02 TEHP=1000. 

47 PDBL.H20 AIB=0.01 FDEL. 02 AIB=0.3 ; 

48 “ 

49 NAHELXSTS 

50 HP='1*B FIT=»0»B P(1)=8 P(2)=7 P(3)=6 P(4)=5 P(5)=4 ; 

51 END 

52 NAHELISTS 

53 TP='1»B C0NDDC='1*B FIT='1»B PBIHT=2 

54 P(1)=8 P(2)=7 P(3)=6 P(4)=5 P(5)=4 P(6)=3 P(7)=2 P(8)=1 P(9)=0.5 

55 T(1)=3000 T(2) =2800 T(3)=2600 T (4) =2400 T(5)=2200 T(6) = 2000 

56 T(7)=1800 T(8)=1600 T(9)=14Q0 T(10) = 1200 T(11)=1000 T(12)=800 

57 ; 

58 END 
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answers for further processing by other procedures. The input data to 
NASAIN is found in file CARDIN. The names, CNDCT, JANAF, FITDAT, and 
CARDIN, are built into the NASA code and, hence, cannot be changed. 

This problem requires two distinct operations, 1) the determination 
of the flame temperatures, and 2) the calculations of the thermodynamic and 
electrical properties. Each of these, in turn, requires three processes: 
1) inputting data, 2) calculating the results, and 3) printing the results. 
Line 31 indicates the simplest GPSAP input for performing these three tasks 
with the NASA code. The MMM statement will call the NASAIN entry, which 
will read the input data on file CARDIN down to the first END keyword. It 
Vjill then call NASA and perform the computations as dictated in the input 
data. Finally, it will call NASAOUT to print a summary of the results, 
which will be the flame temperatures in this ease. 

The second line, 32, in file STRUCIN will perform exactly the same. 
This time, the first call to NASAIN will read down to the next END keyword, 
encompassing the inputs necessary to perform the calculations of the thermo- 
dynamic properties. 

Lines 39 to 58 represent the input data to the NASA code for this 
problem. As stated, the problem can be input directly without any prior 
anlaysis. Lines 39 to 47 represent the FUEL group, giving the coal com- 
position, moisture readjustments, ash rejection corrections, stoichiometry, 
etc. Lines 52 to 57 represent the NAMELIST input. Thus, the HP mode will 
be used with the desired pressures. The END keyword in line 51 indicates 
the end of the data that will be read-in in this call to NASAIN. 

Lines 52 to 56 indicate the NAMELIST input necessary to carry out the 
thermodynamic properties calculations. Line 58 signifies the end for the 
second NASAIN call. 
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Fig. IV-2 Schematic diagram of proposed combination combustor 


B. Example 2 
1 . Problem 

Operate the combination combustor, depicted in Fig. IV-2, of a pro- 
posed experimental MHD facility. The vitiator burns a mixture of No. 2 
fuel oil and preheated oxygen-enriched air. Oxygen is added to the air to 
establish 30 mass % O2. In the combustor, the output of the vitiator is 
combined with additional, room- temperature oxygen, with methane, and with a 
one-to-one mixture of potassium carbonate and water. Potassium is to form 
1 mass % of the total reactant mixture. 

The composition of No. 2 fuel oil will be taken as follows: 

C 86.4% 

H 12.7% 

S 0.7%, and 

0 0 . 2 %, 

with a higher heating value of 10,877.3 cal /kg. 

For a number of different air preheat temperatures and a number of 
different (total) oxygen enrichment levels, it is desirable to know the 
relative amounts of oil and methane necessary to achieve stated flame 
temperatures at 3 atm pressure, 7% system heat loss. Specifically, for 
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T preheat = 700, 800, 900, 1000 K, and 
02 enrichment = 0.30, 0.35, 0.40, 
we wish to obtain, 

T flame = 2900, 3000, 3100 K. 

It is required, further, that these temperatures be achieved with the 
minimum possible amount of oil, i.e., an optimal stoichiometry will be found 
simultaneously. 

2. Analysis 

As will be seen, the availability of G®SAP code-creation and optimiza- 
tion capabilities allows this entire study I:o be performed in a single 
computer job. A small amount of preliminary analysis is necessary, however. 

Note, first, that only the final state of the combustion products is 
of interest. This final state is unaffected by the fact that combustion is 
carried out in a two-stage process. We will, accordingly, merely compute 
the equilibrium state that exists when all reactants are combined. Second, 
note that the thermodynamic properties under investigation are intrinsic 
quantities, unaffected by the total reactant mass. Therefore, all reactants 
will be scaled to 100 grams of methane. 

The heating value for methane is not stated and, presumably, is not 
known. Conversely, the heating value for fuel oil is given, and we wish 
to avoid calculating the corresponding 'Hf. ' As now composed, subroutine 
SETUP (which translates the FUEL groups into NASA input arrays) is nther 
awkward for this sort of FUEL/REACTANTS mixture. Therefore, it is easiest 
to input only the fuel oil through the FUEL group, and to adjust the amounts 
of air, O2 (hot and cold), K2CO3, and H2O explicitly from the input 
structure. 
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The question that must be answered is, given, 
CH = 100 grams, 
oil = X grams, 
stoichiometry ratio = ST, 

O2 enrichment = 02_ENRICH, 

T-preheat = AIRT, 

seed fraction = 0.01 , and 


K2CO3/H2O = 1.0, 

find the weight of KgCO^, H2O, O2 (hot), air, and O2 (cold). Once stated 
the problem is answered fairly easily. The total weight of oxygen neces- 
sary (air and pure oxygen) is 


02T0T = ST * 


ro+ 100 . 0.864*X . 1 0.127*X 

jwrm 12 7011 ■■ 2 2.'0i6 ■ 

+ j 270 ^ ] * 32 - 0.002*xj 


= ST * [398.928 + 3.315 * X], 


where, in trve first entry, the four terms within the square brackets are, 
in order, the moles of methane and the moles of C, H, and S in the oil; the 
last term is the mass of oxygen already found in the fuel oil. From this, 
the weight of the air is 

AIRW = 02T0T/02 ENRICH. 

The hot 02-plus-air is to contain 30 mass % O2; therefore the weight 
of hot O2 is 

02_H0T = 0.3*AIRW, 
and the weight of cold O2 is 
02_C0LD = 02T0T - 02_H0T. 

Finally, the weight of potassium carbonate and of its water of hydra- 
tion are determined by the requirement that potassium, which constitutes 
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0.5658 of the weight of K2CO3, is to comprise 0.01 of the total weight. 
The seed weight and the weight of H2O are, thus, each 

SW = 0.01*(100 + X + AIRW + 02_H0T + 02_C0LD)/0.5658 . 

These manipulations are, of course, slightly tedious. It should be 
noted, however, that this is all that is required for the entire study with 
the present form of the code. Under the old regime also, all the same 
steps were required. In addition, 'Hf would have to be calculated. 
And, finally, all formulae would have to be evaluated for each set of the 
X, ST, and 0^ _ENRICH to be studied. 

3 . Input 

Only the input structure on file STRUCIN and the input data to NASAIN 
on file CARDIN will be explained. The rest of the JCL is similar to that 
for Example 1. Again, the input, shown in Figure IV-3, will be discussed 
on the basis of the line numbers shown on the left. 

The REACTANTS cards, introduced by line 52 and contained on lines 53 
to 58, are the same as those in the previous version of the code (except 
that separation into 'fuels' and 'oxidants' has been eliminated). The 
CH^, K^COg, H2O, preheated O2, air, and cold 0^ are described in turn. 
These will be known to NASA as reactants 1 through 6. The relative weights 
of all but reactant 1 (CH^), as well as the enthalpy and reference tempera- 
ture of reactant 5 and the reference temperature of reactant 4, are dummies, 
because they will be changed by the MAIN Program (these numbers could, in 
fact, be deleted). The blank line (126) is necessary to indicate the end of 
the REACTANTS group. 

The FUEL group, lines 60 through 63, defines the fuel oil and the 
relative weights of carbon, hydrogen, and oxygen. Because STOICH and 
SEED FRAC default (see Section III) to 1.0 and 0.01, respectively, it is 
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Figure IT-3 
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20 
21 
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23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 


//OHE.STBDCIN DD • 

HHH N&SKIH ; 

EBB h ilBT = 700 800 900 1000 ; 

AIBH=-2059.7'.+»IBT* (6.713 + 

*IHT* (2.3485E-4 + 

AIBT* (3.823E-7 - 
AIBT*1.174F-10))); 

BBB B 02 EHPICH = 0.3 0.35 0.4 ; 

BBB C ~T0 = 2900 3000 3100 ; 

BBB D X 100 10 500 

ST 0.9 0.75 1.15 ; 

02T0T=ST*(398.928 + 3.315+X); 
AIBW=02_ENBICH *02T0T; 

O2_HOT=0. 3*ATPW; 

02 COID=HAX (0.001, 02T0T-02 HOT); 

SK= (100 + X+AIPB + 02 H0T+02_C0LD) /54.58; 
AAA 'FOEI.TOT HT • X 


•PFCSTS (2) • 

SB 

•PECBTS (3) • 

SB 

•PFCBTS (4) ' 

02_H0T 

• ETFHP(4) ' 

AIET 

•PFCBTS (5) ' 

AIBB 

•EH'rH(5) • 

AIBH 

•P7’5BP (5) • 

ATBT 

•PECBTS (6) • 

02 COLD 

•HP^ 

•O^B 

• !?• 

• 1 •B 

•FIT-^ 

•O^B 

•CONDOC 

•O^B 

•PPINT^ 

0 

•BEFH^ 

0.0 

•P(1)' 

1.0 

•T(1) • 

300.0 


HH5 NASA ; 
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AAA •HP^ 


37 

1 B 

•O^B 

38 

•P (1) • 

3 .0 

39 

•T(1) • 

0. 93 

40 

• PEFH^ 

NASAP->A(4, NASAP->IFIN) 


41 NHH NASA ; 

42 CCC D NASAP->A (2,1) = TO ; 

43 EEP D X HAXIT=150 ACC=1.0E-2 DEL=1.0E-3 ; 

44 POT SKIP E3IT('AIET=',AIBT,«ONO='OHO,'TO=',TO, 

45 'OIL WT=',X,'ST0ICH=',ST) 

46 (5 (A,E(12,5) ,1(5))); 

47 HHM NASAODT ; 

48 EEE C ; 

49 EEE B ; 

50 EEE A ; 

51 //CABDIH DD * 

52 REACTANTS 


53 

C 1.00 

H 

4.00 

00 

1.5000 

G 

298.1 

F 

54 

K 2.00 

C 

1.00 

0 3.00 

18.902 

S 

298.15 

F 

55 

H 2.00 

0 

1 .00 

CO 

18.902 

L 

298.15 

F 

56 

0 2.00 



00 

238.177 

G 

1073. 5 

0 

57 

N 1.5618 

0 

0.4196 

ABO. 0093 C 0.0003 

693.489 

5730.34 G 

1073.0 

0 

58 

0 2.00 



00 

0.001 

G 

298.15 

0 

59 









60 

FUEL 








61 

FUEL . C= 

= 86 

.4 FUEL.H=12.7 FUEL.S=0.7 

FUEL. 0=0. 2 





62 FOEt.STOICH=0.0 FDEI.SEED PBAC=0.0 FUEL .T=298 . 15 

63 FDEL.HHT=10877.3 ; 

64 

65 NAHELISTS 

66 PBIHT=1 ; 

67 

68 END 
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necessary to zero them explicitly to avoid introducing additional air, 
oxygen, and seed. Input, here, is terminated by the semicolon; the blank 
line (64) is cosmetic only. 

Nearly all of the information input in the NAMELISTS group will be 
controlled by the MAIN Program. An exceedingly ab^'reviated NAMELISTS is 
included to set the defaults and to establish the desired amount of output 
on the 'production' runs. The END card returns control to the MAIN program. 

The MAIN Program, which will perform the desired study, is described 
by the GPSAP input structure, lines 2 to 50. Line 31 calls the input 
module, NASAIN, which reads file CARDIN. Note that the reactants cards are 
read directly into arrays NAME, NUM, PECWTS, ENTH, FAZ, and RTEMP. The 
information from t‘ie FUEL group, on the other hand, is not translated into 
these arrays at that time. 

The three BBB loops-- lines 3 to 50, 8 to 49, and 9 to 48— sweep the 
desired parameter space in air temperature, oxygen enrichment, and flame 
temperature. Lines 4 to 7 calculate the air enthalpy for the specified 
temperature. The optimization loop (lines 10 to 43), carried out for each 
(AIRT,02ENRICH,T0) combination, is actually the heart of the computations. 
The BBB statement in lines 10 and 11 states that an optimizaiton study is 
to be performed and that the independent variables are X (the weight of 
oil) and ST (the overall stoichiometry). The starting value of X is 100 
(one-to-one CH and oil) and will be varied automatically between 10 and 
500; ST begins at 0.9 and remains in the interval 0.75< ST < 1.15. For 
given values of X and ST, the amounts of air, hot oxygen, cold oxygen, 
seed, and water of hydration are calculated in lines 12 to 16. 

The AAA statement, lines 17 to 34, alters a large number of parameters 
within the next calculational entry called— here, the NASA call in line 35. 
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The total weight of the fuel group (the oil) is adjusted to X in line 17. 
The seed weight, SW, is assigned to the potassium carbonate (line 18) and 
to the water (line 19). Reactant 4, the preheated oxygen, is assigned 
weight 02_H0T and temperature AIRT; the air (reactant 5) is given weight 
AIRW, enthalpy AIRH, and temperature AIRT, and the cold oxygen weight, 
02JC0LD, is assign;^'d to reactant 6. Note that a reactant may not have 
negative or zero relative weight; accordingly, 02_C0LD has been adjusted to 
remain > 0.001. 

Lines 26 to 34 define a typical reference enthalpy run. TP is 
true ('I'B), because the enthalpy is to be derived at given pressure 
(1 atm) and temperature (300 K); FIT and CONDUC are false, because poly- 
nomial fitting is not required at this point and conductivity is not 
needed. REFH is zeroed, so that the derived enthalpy is not distorted in 
successive iterations by the reference value derived in the preceding 
iteration. The pressure P(1) is assigned to 1 atm, and the temperature 
T(1 ) to 300 K. 

After the reference run in line 35, the AAA statement in lines 36 to 
40 establishes an HP run at 3 atm and 1 % heat loss. The value of the 
enthalpy at the last calculational point from the previous NASA call (i.e., 
at P = 1 atm, T = 300 K) is assigned to REFH. 

In line 41, NASA is called to calculate the flame temperature at 
3 atm, 7% heat loss. Line 42 states that the optimization in this BBB loop 
will be carried out under the constraint that this flame temperature be 
the given value TO. The EEE statement in line 43 identifies the end of the 
optimization loop and states that X and ST are to be varied so as to 
minimize X subject to the above-stated constraint on flame temperature. 
MAXIT identifies the maximum number of times that lines 10 to 43 may be 
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executed (for each value of AIRT, 02_ENRICH, and TO) in attempting to 
perform the optimization. DEL is the relative step size the optimizer uses 
for computing derivatives, and ACC is the minimum accuracy to which the 
constraint and optimization are to proceed. These last two are decreased 
from their defaults (10“® and 10”^), both to save computation time and to 
avoid numerical difficulties associated with the fact that NASA is a 
single precision routine. 

Line 47 then calls the NASAOUT entry to print the thermodynamic 
properties table, which will consist, in this case, of only one line, repre- 
senting the properties at the last call to NASA. 
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THERMODYNAMIC EQUILIBRIUM PROPERTIES AT ASSIGNED 
TEMPERATURE AND PRESSURE 
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THER^!ODYNAHIC EQUILIBRIUM PROPERTIES AT ASSIGNED 
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THERMODYNAMIC PROPERTIES SUMMARY 
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